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PREFACE 
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by  ARO,  Inc.,  AEDC  Group  (a  Sverdrup  Corporation  Company),  operating  contractor  of 
the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  Project  Number  P32I-86. 
The  project  manager  was  Capt.  K.  H.  Leners,  USAF.  The  manuscript  was  submitted  for 
publication  on  October  22,  1980. 

The  authors  wish  to  express  their  sincere  appreciation  to  H.  N.  Glassman  for  his  work  in 
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1.0  INTRODUCTION 

Several  measurements  of  nitric  oxide  (NO)  concentration  in  the  exhaust  of  jet  engines 
and  combustors  by  probe  sampling  and  ultraviolet  (UV)  resonance  absorption  have  been 
reported  (Refs.  1  through  4).  The  optically  determined  values  of  concentration  have 
appreciably  exceeded  the  probe  determined  values.  In  the  present  reexamination  of  the 
theoretical  basis  for  determining  the  NO  concentration  from  absorption  measurements,  the 
values  previously  reported  have  been  redetermined.  This  report  presents  revisions  to  the 
original  theory,  determination  of  parameters  to  be  used  in  the  theory  from  new  experimental 
data,  and  the  results  of  redetermination  of  the  NO  concentration  for  several  cases. 

The  UV  resonance  absorption  method  for  NO  concentration  measurement  involves 
matching  a  calculated  transmittance  at  the  principal  bandhead  of  the  (0,0)  7-band  of  NO 
through  a  medium  of  known  temperature,  pressure,  and  major  species  concentration 
distribution  to  a  measured  transmittance  of  the  NO  (0,0)  7-band  emitted  from  a  gas 
discharge  tube.  The  calculation  is  made  through  use  of  a  line-by-line  radiative  transfer 
model  (Ref.  5).  Dodge  and  Dusek  (Ref.  6)  recently  reported  both  discrepancies  in  the 
original  model  and  some  less  significant  computer  programming  errors  which  they  found 
while  attempting  to  use  the  computer  code.  These  discoveries  led  to  a  revision  of  the  original 
model  and  redetermination  of  critical  parameters  (Ref.  7).  It  is  necessary  to  reconsider  the 
previously  reported  data  in  light  of  the  revised  model  and  newly  determined  parameters. 

The  radiative  transfer  model  depends,  critically,  on  two  source  parameters  and  two 
parameters  characteristic  of  the  absorbing  molecule.  The  source  parameters  are  the  peak  line 
intensity  distribution  and  the  Doppler  half-width  of  the  lines  emitted  by  the  electric 
discharge  source.  The  absorber  parameters  are  the  absorption  coefficient  at  line  center  and 
the  collisional  broadening  line  width.  The  discrepancy  in  theory  (Ref.  5)  was  contained 
within  the  expression  for  the  absorption  coefficient,  namely,  in  the  treatment  of  the 
rotational  partition  function  and  normalization  of  the  band  oscillator  strength.  The 
reported  broadening  parameter  was  experimentally  determined,  based  on  the  validity  of  that 
model  and,  thus,  was  in  error.  The  source  line  intensity  distribution  was  measured  in  the 
previous  work,  and  the  source  translational  temperature  was  assumed  to  be  slightly  greater 
than  room  temperature. 

In  the  present  work  the  theoretical  model  has  been  revised  to  include  (1)  the  correct 
formulation  of  the  state  population,  (2)  the  correction  of  numerical  errors,  and  (3)  the 
incorporation  of  the  latest  values  of  NO  7-band  line  positions  and  oscillator  strengths  (Ref. 
7).  A  laboratory  program  has  also  been  conducted  that  involved  using  two  absorption 
cells  —  a  high-temperature  cell  and  a  high-pressure  cell  with  respective  path  lengths  of  75.57 
cm  and  51.11  cm.  The  high-temperature  cell  was  operated  over  a  temperature  range  of  from 
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295  K  to  750  K,  and  the  high-pressure  cell  was  operated  at  pressures  ranging  from  0.0  to  3.4 
atm. 

The  source  temperature  was  determined  from  these  data  and  with  the  aid  of  the  corrected 
model.  With  this  new  value  of  source  temperature,  the  model  was  used  to  determine  that 
value  of  the  collisional  broadening  parameter  which  best  fit  the  absorption  data  over  a  range 
of  pressures  and  temperatures.  The  results  of  this  laboratory  study  have  shown  a  nonlinear 
temperature  dependence  of  the  collisional  broadening  not  previously  reported.  Also,  these 
data  prompted  a  reexamination  of  the  relationship  of  pressure  and  temperature  to  the 
broadening  parameter  (fl¬ 
optical  measurements  of  nitric  oxide  concentrations  have  been  previously  reported  for  a 
number  of  combustion  sources.  The  measurements  involving  engine  combustors  and  jet 
engine  exhausts  have  been  recalculated  using  the  corrected  model  and  the  new  parameters; 
the  results  are  included  herein.  A  comparison  of  these  old  and  new  optically  determined 
values  reveals  some  differences,  but  the  original  large  discrepancy  with  probe  measurements 
/  remains  for  high-velocity  combustor  exhausts  (Ref.  2)  and  jet  engine  exhausts  (Refs.  1  and 
4).  It  should  be  noted  that  the  results  of  applying  the  newly  found  temperature  dependence 
of  the  collisional  broadening  to  the  burner  and  combustor  measurements  of  the  interagency 
NO  program  (Refs.  6  and  7)  have  not  been  included. 

2.0  REVIEW  OF  MODEL 


2.1  THEORY 

A  theoretical  model  for  the  transmission  of  Doppler-broadened  resonance  line  radiation 
through  absorbing  media  with  combined  Doppler  and  pressure  broadening  was  given  by 
Davis,  McGregor,  Few,  and  Glassman  (Ref.  5)  with  application  to  the  7-band  of  NO.  The 
formulation  given  in  Ref.  5  was  incorrect  in  that  (1)  Hund’s  case  (b)  statistics  (Ref.  8)  were 
used,  (2)  the  spin  splitting  was  not  properly  taken  into  account,  and  (3)  the  normalization  of 
the  line  oscillator  strength  was  inconsistent  with  the  literature  (Refs.  9  and  10).  The  method 
of  calculation  of  the  partition  function  in  the  model  has  now  been  changed  to  conform  to 
the  angular  momentum  coupling  intermediate  between  Hund’s  cases  (a)  and  (b).  A  brief 
review  of  the  theory  is  presented  in  the  following. 

The  transmission,  r,  of  an  isolated,  Doppler-broadened  emission  line  through  a 
homogeneous  medium  of  length,  (,  is  given  by 


6 


A  E  D  C-T  R  -80-65 


where  I®?  is  the  intensity  of  the  source  line  at  line  center,  pf  is  the  frequency  at  line  center, 
(Asj'j)D  is  the  half-width  of  the  Doppler-broadened  line,  and  k(v)  is  the  monochromatic 
absorption  coefficient  of  the  medium.  The  half-width  of  the  source  line  depends  on  the 
source  temperature  and  is  given  by 


2  in  2  /<Ts 


M  c  2 


(2) 


where  Ms  is  the  molecular  weight  of  the  emitting  molecule,  Ts  is  the  source  temperature,  x  is 
Boltzmann’s  constant,  and  c  is  the  speed  of  light. 


The  absorption  coefficient,  k(p),  contains  contributions,  kj(j>),  from  all  lines  near  the 
frequency,  i'j0.  Each  of  these  lines  will  be  broadened  both  by  the  Doppler  effect  and  by 
collisions.  The  absorption  coefficient  of  the  ith  line  is  given  by  the  Voigt  profile, 


ki  (,,)  =  vVt  f 


(f  e-y ' 


("  i->) 


dv 


(3) 


where  kf  is  the  centerline  Doppler  absorption  coefficient  of  the  ith  line,  and  the  broadening 
parameter,  a,  is  a  parameter  proportional  to  the  ratio  of  the  collision-broadening  half¬ 
width,  (Aaj/j)L,  to  the  Doppler-broadened  half-width,  (A^Od,  in  the  absorbing  medium.  The 
broadening  parameter,  a,  equation  is  as  follows: 


Cf 


(Aa<'  j  )  I) 


The  parameter  w\  in  Eq.  (3)  is  the  Doppler  function 


(4) 


('V',)  n 


^[~in  2 


(5) 


It  can  be  shown  (Refs.  6  and  11)  that 

Ct  =  C  p/T  n 


( 


vv; 


where  p  is  the  pressure,  T  is  the  temperature,  C  is  a  constant  for  any  given  mixture  of 
broadening  species,  and  n  is  the  temperature  dependence  exponent.  The  appropriate 
absorption  coefficient  in  Eq.  (1)  must  include  the  contributions  due  to  all  absorbing  lines 
that  overlap  the  source  line,  j;  thus, 
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Hu)  =  £  k  .(„) 

i 


~  Qe^2 

•  ff2+  (a.;  -y)  2 


(7) 


There  may  be  several  emission  lines  within  a  given  instrument  bandpass,  Aiq.  The 
contributions  of  each  of  these  lines  must  be  summed  to  obtain  the  transmission  within  the 
given  instrument  bandpass.  Thus,  the  transmissivity  within  the  bandpass,  Av\  is  given  by 
substituting  Eq.  (8)  into  Eq.  (1),  summing  over  the  emission  lines,  j,  and  dividing  by  the 
total  intensity  of  the  source  within  the  bandpass  as  follows: 


(8) 


where  k(j>)  is  given  by  Eq.  (8). 


The  line  center  absorption  coefficient  for  a  single  transition  out  of  the  12  possible 
transitions  from  a  given  J ' '  is  given  by 


to  -  Ur24n  i n  2  i  N(v"J"K"p)  f(n’  v’  J1  K1  n"v"J"K") 

1  [  ".c*  7  (VJd  <9> 

where  i  refers  to  the  particular  n'  v'J'K'n"v"J"L"  transition  in  the  A2E ,  x2  II  v-band 
system  of  NO.  The  constants  in  parentheses  are  the  electron  charge,  e,  the  electron  mass, 
Me,  and  the  speed  of  light,  c.  The  population  of  the  v"  vibrational  level  in  the 
J'  'K'  'rotational  sublevel  and  one  of  the  parity  levels,  p,  (+  or  -)  of  the  211^  state  in  the 
intermediate  angular  momentum  coupling  case  between  Hund’s  cases  (a)  and  (b)  is  given  by 


j  -  G(v")  +  Fm(J") 

N(v  M"K"nip) _ (2 J "  +  1)  exp  / _ kT  _ 

^  I  o  I  a  1  2 

2 1  ]  rxp  r  -g(v)-_i  I 

v  l  L.  kT  J  m~ 

where  T  is  the  temperature,  G(v'')  is  the  vibrational  energy  (cm1),  and  Fm(J'')  is  the 
rotational  energy  (cm1)  for  either  the  2n1/2  state  (m  =  1)  or  the  2n3/2  state  (m  =  2).  For 
purposes  of  the  population,  calculation  is  approximated  by  the  Hill  and  Van  Vleck  formulas 
(Ref.  8): 
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K,  =  B, 


(J  +  V2)2  -1-V2  J4O  +  %f  +  y(y  -  4) 


D.,r 


(1  la) 


1  \o  \ 

/  r  \a 

(lib) 

(j  * 

-  1  +  X  y 

]i(j  +  '/2)2  + Y(\  -4) 

where  Y  =  A/Bv  where  A  is  the  splitting  constant,  Bv  is  the  inertial  rotational  constant,  and 
Dv  is  the  centrifugal  rotational  constant.  The  oscillator  strength  for  the  transition  is  given  by 


f(n*  v'  J*  K'  Yi"v"  J"K")  =  jj,K,J"K'i)  J''K")  f(v,  y„)  (12) 

2  J"  +  1  v(v'  v"  ) 

where  (J'K'  J'  'K' ')  and  (v'v' ')  are  the  frequencies  of  the  line  and  band  centers  and  where 
the  normalization  for  the  Honl-London  factors  S(J'K' J'  'K' ')  and  the  normalization  of  the 
band  oscillator  strength  follow  that  of  Tatum  (Ref.  9)  and  are  consistent  with  Pery-Thorne 
and  Banfield  (Ref.  10).  The  Honl-London  factors  are  those  of  Earls  (Ref.  12). 

The  line  positions,  v®,  previously  used  were  those  of  Deeszi  (Ref.  13).  Improved  data 
have  recently  become  available  (Ref.  14),  and  these  data  have  been  correlated  (Ref.  6)  into  a 
consistent  formulation.  These  new  line  positions  have  been  incorporated  into  the  model. 

The  source  parameters  are  the  peak  line  intensity,  I^»  and  the  Doppler  half-width  of  the 
lines,  (As^j)d.  It  is  sufficient  to  determine  only  the  relative  peak  line  intensity  since  any 
absolute  magnitude  will  cancel  in  Eq.  (8).  The  relative  peak  line  intensity  can  be  determined 
by  direct  measurement  for  resolved  lines;  the  intensities  of  the  unresolved  lines  can  be 
determined  from  the  upper  state  population  implied  by  the  measured  line  intensities. 
However,  for  low-pressure,  Doppler-broadened  source  lines  in  the  ultraviolet,  direct 
measurement  of  either  the  half-width  or  the  source  temperature  is  difficult.  Section  2.2 
describes  an  indirect  method  of  determining  source  temperature. 

The  absorber  parameters  are  the  absorption  coefficient  at  line  center,  kf >  and  the 
quantity  A,  related  to  the  coilisional  broadening.  If  the  temperature  and  density  within  the 
absorption  path  are  known,  then  kj)  is  determined  by  Eq.  (10)  if  the  band  oscillator  strength, 
f(v' v' '),  is  also  known.  Thus,  when  kf  is  known,  the  broadening  coefficient,  C,  of  Eq.  (6) 
can  be  determined  from  the  value  of  broadening  parameter,  0,  required  to  obtain  a  match 
between  the  measured  transmissivity  and  the  transmissivity  calculated  from  the  model  based 
on  Eqs.  (8)  through  (12). 
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Successful  application  of  the  model  represented  by  Eqs.  (8)  through  (12)  requires  a 
known  source  line  intensity  distribution,  source  temperature,  band  oscillator  strength,  and 
broadening  coefficient. 

2.2  PARAMETERS 

The  model  represented  by  Eqs.  (8)  through  (12),  together  with  carefully  controlled 
laboratory  experiments,  can  be  used  to  determine  the  source  temperature,  Ts,  and  the 
broadening  parameter,  (3,  when  the  source  line  positions  and  intensities  are  specified  and  the 
band  oscillator  strength  is  known.  Careful  measurements  of  oscillator  strength  for  the  NO 
molecule  have  been  made  with  several  methods  (e.g.,  Refs.  10,  15,  and  16).  Review  of  these 
measurements  shows  that  the  value  [Hassan  et  al.  (Ref.  15)]  of  4.09  ±  0.1  x  1(M  for  the 
oscillator  strength  of  NO  is  the  most  reliable. 

2.2.1  Source  Lines 

The  line  positions  of  the  (0,0)  band  of  the  y  system  of  NO  as  specified  by  Engleman  (Ref. 
14)  have  been  adopted  as  the  most  accurate  available.  The  source  line  intensity  distribution 
of  the  AEDC  capillary  discharge  lamp  has  been  remeasured  with  resolution  sufficient  to 
identify  many  lines  lying  in  each  of  the  independent  branches  of  the  (0,0)  band.  No  changes 
from  the  previously  reported  intensity  distribution  (Ref.  5)  were  found. 

2.2.2  Source  Temperature 

Determining  the  effective  source  temperature  hinges  on  the  fact  that  at  low  pressure  the 
transmissivity  becomes  independent  of  the  other  unknown  parameter,  the  collision¬ 
broadening  parameter.  In  the  low-pressure  limit,  the  combined  Doppler-  and  collision- 
broadened  absorption  profile  [Eq.  (3)]  reduces  to  a  Gaussian  profile  whose  half-width 
depends  only  upon  the  absorber  temperature.  Thus,  if  the  absorber  temperature  and  number 
density  are  known  and  if  the  measured  value  of  oscillator  strength  is  used,  then  the  only 
unknown  in  Eq.  (8)  is  the  Doppler  half-width  of  the  source.  This  half-width  is  directly 
related  to  the  source  temperature  through  Eq.  (2). 

The  experimental  apparatus  described  in  Ref.  5  was  used  to  make  a  carefully  controlled 
set  of  transmission  measurements.  A  sequence  of  transmissivity  measurements  at  the  second 
bandhead  was  made  over  a  range  of  pressures,  with  the  NO  number  density  and  temperature 
held  constant.  In  order  to  determine  the  source  temperature,  the  transmissivity  for  a  given 
number  density  and  temperature  was  determined  from  the  intercept  of  the  transmissivity 
versus  pressure  curve  (Fig.  1)  in  the  limit  of  zero  pressure  (corresponding  to  (3  =  0). 
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Eq.  (8),  this  value  of  transmissivity  was  used  in  the  limit  where  (2  approaches  zero  to 
determine  the  source  temperature.  The  source  temperature  was  adjusted  to  force  a  match 
between  the  calculated  and  extrapolated,  measured  transmissivities.  A  value  of  source 
temperature  of  950  ±  25  K  was  obtained  by  this  method.  Again,  it  should  be  noted  that  the 
source  temperature  is  a  characteristic  of  the  lamp  design  and  is  expected  to  vary  for  different 
designs. 
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2.2.3  Broadening  Coefficient 

In  the  original  work  in  which  the  broadening  parameter,  (f,  was  reported  (Refs.  5  and 
17),  a  simple  collisional  theory  was  assumed  in  which  n  of  Eq.  (6)  was  set  equal  to  unity.  The 
data  scatter,  together  with  the  influence  of  a  mecury  effect  on  the  NO  density  caused  by  the 
presence  of  a  McCleod  gage  connected  to  the  absorption  cell,  obscured  any  nonlinear 
dependence  of  a  on  (1/T).  In  more  recent  work  (Refs.  7  and  18), data  obtained  for  elevated 
temperature  were  insufficient  to  evaluate  the  temperature  dependence;  the  Lorentz 
dependence  (n  =  1)  was  again  assumed.  However,  Dodge  et  al.  (Ref.  6)  had  found  that  their 
transmission  measurements  on  flat-flame,  seeded  H2/O2  burners  at  elevated  temperatures 
could  best  be  correlated  with  a  value  n  =  1.2  in  Eq.  (6).  This  warranted  additional 
investigation  of  the  temperature  dependence  of  the  broadening  parameter. 

A  series  of  experiments  was  performed  in  which  known  concentrations  of  NO  were  used, 
with  N2  as  a  diluent  gas,  over  a  range  of  pressures  (from  0.06  to  3.4  atm)  and  temperatures 
(from  295  to  750  K).  Two  absorption  cells  were  constructed  to  accomplish  this  task.  The 
high-temperature  cell  (quartz)  and  the  high-pressure  cell  (stainless  steel  with  quartz 
windows)  had  optical  path  lengths  of  75.57  and  51.11  cm,  respectively.  In  this  investigation 
many  data  points  were  obtained  from  295  K  room  temperature  to  750  K  for  the 
determination  of  the  temperature  dependence  exponent  n  for  the  NO  7-band.  The  measured 
values  of  transmissivity,  together  with  Eqs.  (8)  through  (12),  were  used  to  determine  the 
value  of  A  for  each  experimental  condition.  The  value  of  A  permits  determining  the 
exponent  n  according  to  A  =  C  P/Tn.  The  value  of  n  is  obtained  by  plotting  fn  A/ p  versus  fn 
T,  where  the  slope  is  equal  to  n  and  the  intercept  is  equal  to  fn  C.  Figure  2  is  a  plot  of 
fn(A/p),  normalized  to  (A/p)0  at  the  fill  conditions  of  the  absorption  cell,  versus  fn  T/T0. 
The  temperature  dependence  found  from  Fig.  2  was  n  =  1.5(±0.2);  i.e.,  the  broadening 
parameter  a  obeys  the  relationship  A  =  C  P/T-15. 

Figure  3  is  a  plot  of  ft  versus  P/T1-5  or  175  data  points  showing  the  least-squares  straight 
line  fit  to  the  data.  The  slope  is  another  way  of  determining  the  value  of  C.  From  the  data 
shown  in  Figs.  2  and  3,  the  value  of  C  was  determined  to  be  3.81  (±0.40)  x  104  K'Vatm. 
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3.0  RESULTS  OF  REDETERMINATION  OF  NO  CONCENTRATION 
FOR  SEVERAL  COMBUSTION  SOURCES 

The  absorption  data  for  the  following  cases  (Refs.  1  through  4)  have  been  subjected  to 
the  revised  theoretical  model  and  computational  procedure:  (1)  YJ-93  Turbojet  Engine  (Ref. 
1),  (2)  Jet  Engine  Combustor  Exhaust  (Ref.  2),  (3)  Turbine  Engine  Combustor  (Ref.  3),  and 
(4)  F101-GE-100  Turbofan  Engine  (Ref.  4).  The  results  are  presented  and  discussed  in  the 
following  paragraphs. 

3.1  YJ-93  TURBOJET  ENGINE  (REF.  1) 

Reference  1  reports  the  results  of  a  measurement  program  whose  purpose  was  to  study 
some  of  the  constituents  of  the  exhaust  plume  of  a  YJ93-GE-3  engine  at  simulated  flight 
conditions.  This  was  the  first  reported  use  of  the  NO  narrow-line  UV  absorption  technique 
in  conjunction  with  a  high-temperature  combustion  process.  The  transmission  data  were 
reduced  to  concentrations  by  extrapolation  of  laboratory  calibration  data  and  the 
assumption  of  Doppler  broadening.  The  line-by-line  transmission  model  (Ref.  5)  had  not 
been  developed  at  that  time.  NO  concentrations  were  measured  at  flight  conditions  from 
Mach  1.4  at  35,000  ft  to  Mach  7.6  at  65,000  ft. 

Table  1  shows  the  results  of  reducing  some  of  these  data  by  obtaining  an  average  static 
temperature  and  static  pressure,  and  using  the  centerline  bandhead  transmission  to  calculate 
an  average  NO  concentration  in  ppmv.  Table  1  gives  the  original  data  (Ref.  1)  and  the  data 
obtained  with  use  of  the  improved  model  and  compares  them  to  an  average  probe  value.  The 
sample  gas  extraction  data  were  obtained  in  the  study  with  a  water-cooled,  constant-area, 
stainless  steel  probe. 

3.2  JET  ENGINE  COMBUSTOR  EXHAUST  (REF.  2) 

Reference  2  considers  NO  emission  from  a  T-56  turbine  engine  combustor  using  JP-4  at  a 
constant  fuel/air  ratio  of  about  0.013.  Pyridene  was  added  to  the  fuel  to  increase  NO 
production  under  certain  conditions.  Table  1  shows  the  results  of  calculating  the  average  NO 
concentration  by  the  method  described  previously.  The  results  of  the  original  model 
(reported  in  Table  1,  AEDC-TR-76-134)  and  of  the  corrected  model  are  compared  to  an 
average  probe  value. 

The  probe  measured  static  temperature  and  static  pressure  profiles.  Transmission 
profiles  were  also  used  directly  in  an  inversion  technique  to  derive  an  NO  concentration 
radial  profile  for  the  two  downstream  stations  reported.  The  details  of  the  original  inversion 
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technique  and  a  revised  technique  are  discussed  in  Refs.  2  and  7,  respectively.  The  original 
and  corrected  data  for  the  stations  3  in.  and  18  in.  downstream  are  shown  in  Figs.  4  and  5, 
respectively.  The  probe  used  was  an  uncooled,  constant-area,  stainless  steel  probe. 

3.3  TURBINE  ENGINE  COMBUSTOR  (REF.  3) 

This  report  documents  the  measurement  of  exhaust  product  species  on  a  14.0-cm-diam 
AVCO  Lycoming  development  turbine  engine  combustor  that  was  loaned  to  AEDC.  Data 
were  taken  at  various  fuel/air  ratios  (0.01  to  0.05)  at  a  station  1 .27  cm  downstream  from  the 
combustor  nozzle  exit.  Table  1  shows  the  results  of  calculating  the  average  NO 
concentration.  As  before,  the  results  of  the  original  model  (Table  5,  AEDC-TR-76-180)  and 
the  corrected  model  are  compared  to  an  average  probe  value.  A  water-cooled,  constant-area 
probe  was  used  to  obtain  the  gas  sample  data  reported. 

3.4  FlOl-GE-lOO  TURBOFAN  ENGINE  (REF.  4) 

This  report  describes  the  measurement  of  pollutant  species  in  the  exhaust  gases  of  an 
F101-GE-100  jet  engine.  NO  and  NOx  absorption  data  were  taken  at  the  intermediate  and 
maximum  afterburn  power  settings.  The  optical  and  probe  measurements  were  made  at  a 
station  10.2  cm  downstream  of  the  nozzle  exit  of  two  engines  of  different  serial  numbers  that 
had  very  similar  performance  characteristics. 

Gas  sample  probe  data  were  obtained  at  several  simulated  flight  conditions  at  only  the 
intermediate  engine  power  setting.  Comparison  between  probe  and  optical  results  can  only 
be  made,  then,  at  this  power  level.  The  optical  and  probe  data  were  not  taken  at  the  same 
altitude,  but  a  probe  profile  was  interpolated  from  the  values  corresponding  to  altitudes 
which  bracketed  the  optical  data  acquisition  point. 

The  inversion  technique  was  reapplied  to  the  data  reported  in  Ref.  4  at  the  intermediate 
power  setting.  Figure  6  compares  the  NO  and  NOx  probe  profiles  to  the  results  of  the 
original  inversion  technique  as  reported  in  Fig.  11  of  Ref.  4  and  to  the  results  of  the 
inversion  technique  using  the  improved  model.  An  orifice  probe  similar  to  the  “proper 
design”  probes  described  in  Refs.  19  and  20  was  used  to  obtain  the  gas  sample  data  reported. 
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Figure  4.  Comparison  of  the  measured  probe  profiles  with  the  results  of  the 
original  inversion  technique  (Fig.  10a  of  AEDC-TR-76-134)  at  a 
station  7.62  cm  downstream  from  combustor  exhaust  duct  exit. 
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Optical  NO  Data  from  AEDC-TR-76-134 
(Combustor  Rig) 


Radial  Distance,  r/R^.  (R^  =  15.2  cm) 


Figure  5.  Comparison  of  the  probe  profile  with  results  of  the  original 

inversion  technique  (AEDC-TR-76-134)  at  a  station  45.72  cm 
downstream  of  combustor  exhaust  duct  exit. 


Table  1.  Example  Results  for  NO  Measurements  Made  Previously 


Condition 

Average 

Static 

Temperature 

K 

Average 

Static 

Pressure, 

psia 

Path 

Length, 

cm 

Measured 
Transmissivity , 
Optical  Method 

Average  NO 
Cone.,  Optical, 
Originally 

Pub.,  (ppm) 

Average  NO 
Cone.,  Optical, 
Improved 

Model ,  (ppm) 

Ratio  of 
Opt^  to 
Probe 

Average  N0X 
Probe 
Reported 

AEDC-TR-79-132 

YJ93- 

-GE-3  Jet  Engine 

Alt  =  65,000  ft 
Mach  =  2.6 

Military  Power 

506 

0.99 

112.0 

0.626 

323 

436 

4.36 

100 

Alt  *  65,000  ft 
Mach  -  2.6 

Max.  A/B 

868 

1.52 

127.0 

0.594 

617 

686 

5.28 

130 

AEDC-TR-76-134 

WPAFB-T56  Combustor 

7.62  cm 

Dnst .  N.E.* 

847 

14.7 

30.4 

0.884 

170 

119 

1.09 

109 

7.62  cm 

Dnst.  N.E. 
(Pyridene) 

847 

14.7 

30.4 

0.875 

190 

129 

1.03 

125 

7.62  cm 

Dnst.  N.E. 
(pyridene) 

847 

14.7 

30-4 

0.830 

270 

181 

i  , 

1.17 

154 

AEDC-TR-76-180 

Avco 

-Lycoming  Combustor  All  Data  Obtained  1.27  cm  Downstream  of  Nozzle  Exit 

Fuel-to-Air  Ratio: 

0.01 

789 

20.2 

6.0 

0.979 

175 

82 

2.73 

30 

0.02 

1,022 

15.0 

6.0 

0.952 

605 

305 

3.0 

102 

0.03 

1,332 

17.1 

6.0 

0.931 

860 

617 

3.25 

190 

0.04 

1,550 

16.7 

6.0 

0.920 

1,220 

953 

4.43 

215 

0.05 

1,724 

18.1 

6.0 

0.910 

1,600 

1,236 

4.68 

264 

*Downstream  of  Nozzle  Exit 
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NO  and  NO  Concentration,  ppm 


r/R^  (R^  =  Nozzle  Radius) 


Figure  6.  Comparison  of  the  reevaluated  NO  and  NOx  probe  profiles 
with  results  of  the  original  and  improved  model 
(AEDC-TR-77-75)  10.8  cm  downstream  of  engine  nozzle  exit. 
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4.0  DISCUSSION 

4.1  PREVIOUSLY  REPORTED  DATA 

The  optically  determined  NO  concentration  data  reported  in  the  literature  have  been 
corrected;  the  more  significant  data  are  presented  herein.  The  data  presented  here  are 
intended  to  indicate  the  magnitude  of  error  in  the  previously  reported  results. 

1.  For  the  YJ93-GE-3  engine  (Ref.  1)  two  data  points  give  a  higher  NO 
concentration  than  the  original  results,  with  a  correction  for  each  of  about  9  and 
46  percent,  respectively.  The  ratio  of  corrected  optically  determined  values  to 
probe  sampling  values  ranges  from  4.4  to  5.3,  compared  to  previously  reported 
values  of  3.2  to  4.7.  It  is  important  to  note  that  these  data  are  averaged  along  the 
path  and  that  the  original  data  were  reduced  by  extrapolation  of  laboratory 
measurements  rather  than  by  correlation  with  the  radiative  transfer  computer 
model. 

2.  For  the  T-56  combustor  (Ref.  2),  all  newly  determined  average  NO 
concentrations  are  smaller  than  those  previously  reported,  with  an  average 
reduction  of  about  24  percent.  The  ratio  of  corrected  optically  determined  to 
probe-determined  average  values  ranges  from  1.03  to  1.17,  probably  within  the 
data  reproducibility.  For  the  data  reduced  by  the  inversion  procedure,  the 
corrected  NO  concentrations  at  centerline  are  lower  than  those  originally 
reported.  The  corrections  are  46  percent  at  the  3-in.  station  and  40  percent  at  the 
18-in.  station.  The  optically  determined  values  now  agree  very  well  with  the 
probe-determined  values. 

3.  For  the  Avco-Lycomrng  combustor,  the  percent  correction  in  the  optically 
determined  values  ranges  from  22  to  31  percent  for  the  five  different  fuel-to-air 
ratios  used.  The  ratio  of  corrected  optical  to  probe  concentrations  ranges  from 
2.7  to  4.6,  whereas  the  range  of  the  originally  determined  ratios  was  from  5.8 
to  6.0. 

4.  Using  the  inversion  technique,  the  percent  correction  for  the  centerline  value  of 
NO  is  27  percent  for  the  F10i-GE-100.  The  ratio  of  corrected  optical  to  probe 
NOx  centerline  concentration  is  1 .52  compared  to  2. 1  as  previously  reported. The 
ratio  of  corrected  optical  to  probe  NO  centerline  concentration  is  2.00. 
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It  is  to  be  noted  that  the  optically  determined  NO  concentration  is  generally  in  close 
agreement  with  the  probe-determined  NOx  concentration  for  the  low-velocity,  low  fuel-to- 
air  ratio  combustion  streams  encountered  from  the  T-56  combustor  rig  in  which  an  orifice 
probe  was  used.  But,  for  the  higher-velocity,  higher-temperature  reactive  plume  of  the 
Avco-Lycoming  combustor,  the  optically  determined  NO  concentrations  are  a  factor  of  2.4 
to  4.7  greater  than  the  probe  value  where  the  data  are  taken  only  0.5  in.  downstream  from 
the  nozzle  exit  (before  any  possible  mixing  can  occur)  and  where  a  constant-area  probe  was 
used.  Notice  also  that  the  ratio  of  optical-to-probe  values  is  highest  at  the  highest  fuel-to-air 
ratios  and  the  respective  higher  temperatures.  The  jet  engine  data  (NO  to  NOx),  also  under 
high-velocity  and  high-temperature  conditions,  have  relatively  high  optical-to-probe  ratios. 
For  the  YJ-93  turbojet  engine,  the  ratio  of  the  average  concentration  values  was  about  4  for 
military  power  operation;  for  the  F101-GE-100  turbofan  engine  operating  at  military  power, 
the  centerline  ratio  was  about  1 .5  for  probe  NOx,  but  in  the  cooler,  nonreacting  fan  region, 
the  optical  and  probe  values  agreed  quite  closely.  Thus,  it  appears  that  plume  mixing  and 
temperature  may  be  extremely  important  factors  in  determining  the  relative  differences 
between  the  results  obtained  with  the  probe  and  the  optical  technique.  Trends  in  the  probe- 
measured  NO  and  NOx  concentrations,  however,  lead  one  to  conclude  that  some  chemical 
effects  take  place  within  the  probe  and/or  sample  line.  Different  probe  designs  gave 
significantly  different  NO  concentrations;  therefore,  probes  must  be  specially  designed  for 
each  application  (Ref.  20). 

5.0  CONCLUSIONS 

Although  errors  were  found  in  the  original  model  for  transmission  of  NO  7-band 
resonance  bands  through  absorbing  media,  they  were  partly  self-compensating  because  of 
the  use  of  the  model  in  determining  the  line-broadening  parameter,  which  plays  a  major  role 
in  the  numerical  calculations.  Thus,  previously  published  values  of  NO  concentration 
determined  by  the  optical  method  in  high-velocity  and  high-temperature  streams  did  not 
differ  greatly  from  values  obtained  with  the  improved  model.  This  was  demonstrated  by  all 
the  cases  (both  high-velocity  and  high-temperature)  examined,  and  the  previously  reported 
large  differences  in  optical  and  probe  sampling  determinations  of  concentrations  in  high- 
velocity,  high-temperature  flows  were  sustained.  For  the  case  of  low-velocity  nonreacting 
flows,  the  probe  and  optical  results  were  very  nearly  the  same.  Also,  in  the  mixed,  fan  flow 
region  of  a  turbofan  engine  exhaust,  the  results  for  the  probe  and  optical  techniques  agreed. 
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NOTE  OF  CORRECTION:  BROADENING  OF  NO  y-BAND  LINESt 

W.  K.  McGregor,*  J.  D.  Few.*  D.  R.  Keefer.S  H.  S.  Lowry.  III.  t  and  M.  G.  Davis§ 


Sverdrup/ ARO,  Iruu,  AEDC  Division,  Arnold  Air  Force  Slation.  TN  37389.  U.S.A. 

(Rttdrtd  30  August  1979) 

Abstract— Errors  in  Ref,  1  ore  corrected.  Using  the  new  data,  the  collisions!  broadening  coefficient  of  NO 
rtund  lines  by  N,  is  2390  (±220)  atnf'-K  (old  value  =  1270=200)  and  the  optical  ooHisioodiameter  a 
13JsO.fi  A  (instead  of  3Js0.3  A).  Values  for  argon  (2l50±2MatnT  -K)  and  COj  (2040s 200 atm  -K) 
yield  the  optical  collision  diameters  of  13.2  and  13.1  A,  respectively. 


INTRODUCTION 

Some  of  us  previously'  determined  broadening  parameters  in  electronic  transition  bands  of 
diatomic  molecules.  Dodge  and  Duselr  have  found  that  our  expressions  for  the  ground  state 
partition  function  for  NO  and  the  normalization  of  the  oscillator  strength  values  were  in¬ 
consistent  with  some  published  work3  and  also  identified  minor  errors  in  the  computer  program. 
We  have  reexamined  and  corrected  the  theory,  made  more  precise  laboratory  measurements, 
and  redetermined  the  broadening  parameter  for  the  y-bands  of  NO. 


REVISION  OF  MODEL 

The  X:n  state  of  the  NO  molecule  belongs  to  the  coupling  case  intermediate  between 
Hund’s  cases  (a)  and  lb)4  for  which  the  rotational  energy  levels  may  be  calculated  ap¬ 
proximately  from  Ref.  5  while  the  rotational  line  strength  (Honl-London)  factors  follow  the 
Earls  formula*  normalized  according  to  Tatum.’  In  the  original  paper1,  Hurd’s  case  (b)  statistics 
were  erroneously  applied.  The  population  density  of  quantum  level  z\  /*,  spin  sub-level  c.  and 
parity  p  is  for  the  intermediate  coupling  case 

,J  -  J-  )  - ,  (I, 

2S{exp[-G(»")/ic71  S  2  (2J*+ l)exp[-F((.T)/*71} 

f-l 

where  G(iO  is  the  vibrational  energy,  T  is  the  temperature,  *  is  Boltzmann’s  constant,  and 
FeU”)  is  the  rotational  energy  given  approximately  by  the  Hill  and  van  Vleck  formulae. 
Equation  (1)  yields  values  of  approximately  half  as  large  as  Eq.  (8)  of  Ref.  I. 

Equation  (9)  of  Ref.  1  does  not  comply  with  the  normalization  in  Ref.  3  which  was  used  for 
the  band  oscillator  strength  the  correct  expression  for  the  line  oscillator  strength  Atkvtk- 
is 
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(The  research  reported  herein  was  performed  at  the  Arnold  Engineering  Development  Cem  (AEDC),  Air  Force 
Systems  Command  by  ARO.  Inc.,  a  Sverdrup  Corporation  Company,  operating  contractor  of  AEDC.  Further  reproduction 
is  authorized  to  satisfy  needs  of  the  U.S.  Government. 
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where  i 'rtcru-  and  *,>•  are  the  frequencies  of  the  line  and  band  centers,  respectively,  and 
Sj'ktk"  is  the  Honl- London  factor.  We  prefer  the  most  recently  reported*  value  4.09  ±  0. 1  x 
10  for  the  band  oscillator  strength  of  the  0,0  band  of  the  NO  y-system  in  place  ofJ 
3.64  ±0.05  x  10r*. 

When  Eqs.  (1)  and  (2)  are  used,  Eq.  (7)  of  Ref.  I  becomes 
r  •  • 

k  2e2V(g,  ln2)  N(i;*/*cp)/'„v  tr-r-r.^ 

**  (££  •  (3) 

where  e  and  m  are  the  electronic  charge  and  mass,  c  is  the  speed  of  light  and  i$  the 
Doppler  line  width.  The  computer  model  was  reconstructed  using  Eqs.  (1H3).  Also,  the 
improved  line  locations  of  Engleman,  Rouse,  and  Baramonte*  were  used  rather  than  those  of 
Deezi,  as  suggested  by  Dodge  and  Duse  k. 2 


EXPERIMENTAL  DATA 

Many  experimental  data  were  obtained  using  the  apparatus  of  Ref.  1.  We  now  believe  that 
some  of  the  NO  densities  reported  in  Ref.  1  were  in  error  because  of  interactions  of  NO  with 
the  Hg  in  the  McLeod  gauge. 

The  effective  temperature  of  the  source  was  measured  in  the  current  work  rather  than 
assumed.  In  the  limit  of  zero  partial  pressure  of  broadening  gas  and  with  a  small  absorber 
pressure  (— 10'  atm),  the  absorber  line  shape  reduces  to  a  Doppler  profile,  and  the  trans¬ 
missivity  of  a  narrow  line  source  depends  only  upon  the  source  line  width.  The  measured 
transmissivity  was  extrapolated  to  zero  broadening  gas  pressure  (see  Fig.  1).  The  source 
temperature  in  the  model  was  adjusted  until  the  transmissivity  at  zero  broadening  gas  pressure 
agreed  with  the  extrapolated  measured  value.  The  result  was  950  ( ±  25)°K. 

Values  of  the  Voigt  line  profile  parameter  o'  required  to  match  measured  and  calculated 
transmissivities  at  the  2nd  bandhead  of  the  0, 0  y-band  were  plotted  vs  PIT,  and  the  slope  was 
determined  for  either  Nj,  Ar.  or  COj  as  the  foreign  gas.  Results  for  Nz  are  shown  on  Fig.  2  and 


Fig.  I.  Least-squares  At  to  die  traramisiivity  at  the  first  bandhead  of  the  0, 0  NO  y-baod  vs  PIT  in  the  limit 

u  the  pressure  goes  to  zero. 
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Eg.  1  Valuer  of  the  spectral  broadening  parameter  (s')  aa  a  function  of  pIT  for  the  (0. 0)  ybaod  of  NO. 


broadening  coefficients  obtained  are  given  as  follows: 

a'=C  PIT 

for  N2i  a' a  (2390  ±  220)  PIT ; 
for  Ar,  a'“(2150±  200)  PIT; 

1  for  COjt  a!  -  (2090  ±  200)  PIT;  (4) 


here  the  value  of  C  has  the  units  atnT’-K. 
For  bimolecular  broadening  collisions,1  ‘ 


<3) 


where  <r*  is  the  collisional  cross  section  for  broadening  of  the  absorbing  molecules  a  by 
molecules  of  species  b,  A  is  the  wavelength,  Pa  is  atmospheric  pressure  at  standard  conditions, 
ir  is  Boltzmann's  constant,  Af«  is  the  mass  of  the  primary  molecule,  and  Mb  is  the  mass  of  the 
foreign  broadening  molecule.  Thus,  for  the  NO  y-band.  we  obtain  the  data  listed  below. 


Broadening  gas 

Cross  section  <I0"U  cm2) 

Optical 

collision  diameter  (A) 

Nj 

556*53 

13.3=0.6 

Ar 

547  *  53 

13.2=0.6 

CO, 

540-53 

13.1=0.6 
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In  our  original  paper,  the  exponent  of  the  constant  in  Eq.  (11)  was  incorrect;  it  should  have 
been  10"  and  not  10".  The  result  was  an  error  in  the  reported  optical  collisional  diameter;  it 
should  have  been  II  A  and  not  3.5  A.  The  values  of  the  collisional  diameter  (4*,  defined  by 
“  ***)  reported  here  are  about  13  A  and  are  in  substantial  disagreement  with  that  of 
Thorsen  and  Badger.11 
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